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Abstract 



We summarize a discovery potential for supersymmetric particles at the 
pp collider of Tevatron with center-of-mass energy y/s = 2 TeV and inte- 
grated luminosity J C dt = 15-30 fb _1 . Any direct search is kinematically 
limited to below 450 GeV/c 2 . We, however, have a unique opportunity to 
test various supersymmetric scenarios by a measurement of the branch- 
ing ratio for the rare decay mode B s — * fi + fj," . Using the background 
estimate in the CDF analysis of B s — > fJ, + (J,~ in Run I, we investigate 
the prospects for studying this mode in Run II. CDF would be sensitive 
to this decay for a branching ratio > 1.2 x 10 -8 with 15 fb _1 (or, if a 
similar analysis holds for D0, > 6.5 x 10~ 9 for the combined data). For 
tan (3 > 30, the B s — > \x + \x~ search can probe the SUSY parameter space 
that cannot be probed by direct production of SUSY particles at Run II. 
An observation of B s — ► with a large branching ratio > 7(14) x 10 -8 

(feasible with only 2 fb _1 ) would be sufficient to exclude the mSUGRA 
model for tan (3 < 50(55) including other experimental constraints. For 
some models, the branching ratio can be large enough to be detected even 
for small tan/3 and large m^ 2 . 



*Representing the CDF and D0 Collaborations. Plenary talk at 10th International Confer- 
ence on Supersymmetry and Unification of Fundamental Interactions (SUSY02), June 17-23, 
2002, DESY, Hamburg, Germany. 



I. INTRODUCTION 



The Fermilab Tevatron collider will define the high energy frontier of particle physics 
while CERN's Large Hadron Collider is being built. The first stage (Run Ha) of the 
Tevatron collider Run II (March 2001-2005) will deliver at least 2 fb* 1 of integrated 
luminosity per experiment at \/s = 2 TeV. Major upgrades of the CDF and D0 de- 
tectors [1,2] have been completed and we have beeing taking data. In the second stage 
(Run lib) of Run II (starting in 2006), we expect 15 fb _1 of integrated luminosity per 
experiment at 2 TeV. Further upgrades of the two detectors are been undertaken. 

Among various features, the detectors have the ability to trigger on displaced ver- 
tices from bottom and charm decays using a precise microvertex detector to enchance 
the Higgs search and the physics with top quarks. Searches for supersymmetry (SUSY) 
are among the main priorities along with Higgs and top physics for Run II. 

Supersymmetry uniquely opens the possibility to directly connect the Standard 
Model (SM) with an ultimate unification of the fundamental interactions. With the 
results on electrweak and strong gauge couplings from CERN's e + e~ collider LEP 
experiments and the top quatk mass at CDF and D0 at the Tevatron, the models of 
SUSY have become more predictive and require a spectrum of new particles below a 
few TeV/c 2 . Thus SUSY represents a natural candidate for the new physics expected 
to occur in the TeV energy domain. 

In this paper, we summarize the prospects of SUSY searches/discovery either di- 
rectly through collider processes or indirectly through rare processes at the Tevatron. 

II. SUSY MODELS 

One of the difficulties in determining predictions of generic Minimal Supersymmet- 
ric Standard Model (MSSM) lies in the large number of new parameters (over 100 free 
parameters) the theory implies. One may consider a theoretical framework to reduce 
the number of free parameters. 

Fortunately SUSY models apply to a large number of different accelerator and 
cosmological phenomena, and a great deal of effort has been involved in recent years to 
use the data to limit the parameter space. Part of the difficulty in doing this resides in 
the success of the model in not disturbing the excellent agreement of the precision tests 
of the SM [3] due to the SUSY decoupling theorems which suppress SUSY contributions 
at low energies. Historically, the absence of flavor changing neutral currents at the tree 
level played an important role in the construction of the SM. They represent therefore 
an important class of phenomena that might show the presence of new physics, since 
the SM and the SUSY contributions contribute first at the loop level with comparable 
size. Thus the decay b — > 57 has been a powerful tool in limiting the SUSY parameter 
space. 

Extensive Monte Carlo (MC) studies were carried out during 1998 on the following 
four topics to maximize the direct SUSY/Higgs searches in Run II: (i) supergrav- 
ity (SUGRA) [4], (ii) Gauge-mediated SUSY Breaking (GMSB) [5], (iii) beyond the 
MSSM [6], and (iv) Higgs [7]. The readers can refer to Ref. [8] for a summary of Run I 



SUSY searches. Among the experimental aspects, we conclude that it is important 
to have excellent triggering/tagging and identification for 6's, r's, 7's as well as e's 
and /i's. Thus we develop (A) low px lepton+track trigger, where the "track" object 
can be electron, muon, or hadronically decaying r-lepton, (Bl) >2-jets trigger + frp 
(> 25 GeV) and (B2) >2-jets trigger + (> 20 GeV with b/c tagging), instead of 
inclusive pr trigger, (C) better trigger/identification for prompt /displaced photons. 
Trigger Bl is extremely useful for reliable parametrization of high frp distribution due 
to QCD events to reduce the systematic uncertainty [4]. 

We first consider mininal SUGRA (mSUGRA) and GMSB frameworks as examples 
of direct searches for SUSY production, which characterize the experimental triggers 
and analyses in 6's, r's, 7's. Details of all SUSY models and experimental prospects 
can be found in Refs. [4-6]. We then take the decay B s — > fi + fi~ as an example of 
powerful tool in limiting SUSY parameter space at the Tevatron. 

III. DIRECT SEARCHES 

A. Testing mSUGRA 

The mSUGRA model [9,10] depends on four parameters: m (the universal scalar 
mass at Mq), ^1/2 (the universal gaugino mass at Mq), A (the universal cubic soft 
breaking mass at Mq), and tan/5 (the ratio of two SUSY Higgs vacumm expectation 
values at the electroweak scale). In addition, the sign of \i (the Higgs mixing parameter) 
is arbitrary. With R parity invariance, the lightest neutralino (x?) is assumed to be 
the lightest suppersymmetric particle (LSP) and it is bino-like and stable. The Xi then 
would pass through the detector without interaction. We fix A = and the sign of \i 
to be positive for simplicity, otherwise stated. Here the ISAJET sign convention for /z 
is used. 

1. Chargino- Neutralino Associated Production 

The trilepton signal arises when both the lightest chargino (xf) an d the next-to- 
lightest neutralino (xV) decay leptonically in pp — > xfx® + X. An initial study of the 
final state of trilepton (£££) plus frp for high luminosity (>10 fb _1 ) at the Tevatron 
was made in Ref. [11] for direct X\X\ production. Here tit = eee, ee/z, e/z/i, /i/x/i by 
excluding electron or muon from tau (r) lepton leptonic decay. Further studies were 
also made for pp — > 3£ + X, including all SUSY production processes (e.g., £v) and 
its decays [12-14,4,15]. Electron or muon from the leptonic decay mode of r lepton 
were accepted. This requires CDF and D0 experiments to trigger and identify the 
leptons with low p T (>10 GeV/c). They also improved the SM background calculations 
including effects from W*, Z*, and 7*. 

The studies in Refs. [12,14,4,15] have included hadronically decaying r lepton as 
well. This requires both experiments to trigger and identify the r leptons with low 
Pt (^10 GeV/c). With high luminosity in Run II, an inclusive lepton trigger becomes 
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FIG. 1. Regions of the m^j-mo plane where the trilepton events should be detectable at a 
level of 5a significance for tan/3 = 5 (left) and tan/3 = 50 (right) [15]. Three areas are shown 
for an integrated luminosity of 30 fb _1 (magenta), 10 fb^ 1 (blue) and 2 fb _1 (green) and from 
top to bottom, respectively. The large red regions are excluded by theory and experiment. 
Dashed lines represent the SUSY contribution to the muon anomalous magnetic moment 
(in units of 10 -10 ) and the dotted lines are iso-mass contours of the lightest neutral Higgs 
boson. The solid contour (only for tan/3 = 50) indicates the prediction for the branching 
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more difficult because of such a large trigger rate. Thus a generic dilepton trigger, 
namely lepton + track trigger, is neccessary [4]. 

In those studies, the final states of 3£ + fir, £ £ + th + fir are found to be 
the best channels for the study of chargino-neutralino associated production. Trig- 
ger A, mentioned earlier, will play a key role to maximize the experimental sensi- 
tivity in this channels. Figure 1 shows 5a discovery reach in the trilepton channel 
in mSUGRA for small tan/3 = 5 and large tan/3 = 50. We will be sensitive upto 
mi/a ~ 250 GeV/c 2 if m < 200 GeV/c 2 and low tan/3 (e.g., 5). With large tan/3 
(e.g., 50), upto m 1/2 ~ 200 GeV/c 2 if m > 500 GeV/c 2 . It should be noted that the 
CDF and D0 analyses [16,17] of the trilepton channel in Run I (about 100 pb _1 ) were 
limited to tan/3 = 2, /x < within the MSSM framework. The chargino mass upto 
about 80 GeV/c 2 (m^ = 75 GeV/c 2 ) was excluded for tan/3 = 2, /j < and A = 
if Mq = Mg within mSUGRA [18]. 

2. Gluinos and Squarks Production 



Gluinos (g) and squarks (q) are pair-produced at the Tevatron. Within the 

Thus, there appear two representative parame- 



mSUGRA framework, Mg > 0.85 Mg 
ter regions in terms of the production: (i) 



gg, gq, and qq productions where Mg 



M~ g 



and (ii) gg production where Mg ^> Mg. 

In most of the parameter space accessible at the Tevatron, the left-chiral squark 
dominantly decays into a quark and either a xt or a X.2- Those branching ratios are 
B(q L - q'xt) - 65% and B(q L - qx° 2 ) ^ 30%. Since M~ - M~± > M~± - M~ , 

^ A.1 A.1 Al 

the jet in the g L — > g'%f (or g L — > decay likely has larger than those in 

the xf - >■ 9<fXi (o r X2 qq'Xi) decay [19]. Similarly, at least one jet in the gluino 
decay (g — > qq'xt or 99X2 through a real of virtual squark) has large E T . Thus, pair- 
produced squarks and gluinos have at least two large- jets associated with large fir- 
Furthermore, the jet multiplicity tends to be larger for events with gluino than with 
squark. The final state with lepton(s) is possible due to leptonic decays of the xt 
and/or The branching ratio to the final state with two or more leptons strongly 
depends on the value of tan f3. This leads us to look for SUSY events with final states 
of from jets + fir and \t + jets + |r = e or /i) [20]. We restrict the parameter 
space so that lighter third generation squarks (b\ and t\) remain heavier than the xt 
and the x°- 

In the jets + fir channel, for example, an optimization of cuts could be made on Nj, 
fir, and M S2 (= fir +E T 1 +E 3 T 2 ) [20]. The final selection cuts are: (a) Nj > 4; (b) veto 
on isolated leptons (e or fi) with p T > 15 GeV/c; (c) fir > 100 GeV; (d) Acj). * > 30°; 

(e) Ms 2 > 350 GeV. The SM background sizes are estimated to be 25 fb for it events, 
38 fb for W/Z + jets events, 1 fb for diboson process, and 9 fb for QCD events, totaling 
73 fb. Figure 2 is the significance as a function of Mg where Mg ~ Mg at tan /3 — 3, 
10, and 30 (/i > and A — 0). The strongest reach in 5a significance is 410 GeV/c 2 
(m 1/2 ~ 160 GeV/c 2 ) for 15 fb" 1 . There is no significant tan (3 dependence. This can 
be compared to 280 GeV/c 2 for 100 pb~\ 360 GeV/c 2 for 2 fb" 1 , and 440 GeV/c 2 for 
30 fb" 1 . This channel has huge QCD background events, so that Trigger Bl will play 
a key role to minimize systematic uncerainty in understanding the size of the QCD 
events. 

In the 1£ + jets + fir channel, the gluino mass limits are less stringent than those in 
the jets + fir channel [20], except for smaller mo values. It will be essential to combine 
limits from the two channels to maximize the sensitivity in Run II. 

It should be noted that the D0 and CDF analyses of the jets + fir channel within 
the mSUGRA framework (A = 0) in Run I were limited to tan/3 = 3, /i < [21,22]. 
The stringent lower limit on the gluino mass at 95% C.L. is 300 GeV/c 2 (mi/ 2 ^ 
130 GeV/c 2 ) for tan/3 = 3 and /i < if Mg = Mg [22]. The D0 analysis of the le + 
jets + fir channel in Run I was also limited to tan/3 = 3, fi < [23] and the gluino 
mass limit was less stringent than those in the jets + fir channel. 

3. Stop and Sbottom Production 

A large mixing angle 9^ between the superpartners of the left-chiral and the right- 

chiral stop quarks, i L and Ir respectively, form two squark mass eigenstates where 
M- tl < Mi 2 . The ti could substatially be lighter than other squarks. mSUGRA with 
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FIG. 2. Significance as a function of Mg (Mz — Mg) f° r tan/3 = 3 (filled circles), 10 
(down triangles), and 30 (open circles) in jets + tfsr channel (15 fb _1 ) at the Tevatron [20]. 



smaller tuq and/or larger \Aq\ would have a light stop in their spectrum. In contrast, 
lighter sbottom (pi) can appear only at small m and small m 1( / 2 in the mSUGRA 
models. In addition, they are always accompanied by light stops, except tan/3 > 20 
and /j, < 0. 

For general studies, we simply assume eithe b\ or t\ is the lightest squark. A 
comprehensive study on prospects of those searches can be found in Refs. [24,4] . Decays 
studied for the i\ or b\ in Run II are: (i) b\ —> bxi, (ii) h — > cxI, (hi) h — > bW(b£i/), 
(iv) t a - bxf, followed by xt - ^ W Xi -> ^X?, (v) t~i ^ 6VFx?. 

Table I is a summary of the maximum sensitivity in the searches. In Run II, we 

+ (*) 

should also consider the case where the braching ratio for t\ — > bxf —> brv is nearly 
100%. Trigger A should enhance the sensitivity of the serach. It should be noted that 
the above studies are based on bino-like LSP. For higgsino-like LSP, the search strategy 
needs to be modified and its prospects can be found in Ref. [24]. 

The LEP limits on xt an d i masses [25] leave the decays of (a) b\ — > bxi [26,27], 
(b) tj -> cxi [26], (c) ?! -> biv (three-body decay, M 9 = M w ) -> Uvx\ [28], and 
(d) t\ —>■ blvx\ (four-body decay, My = Mr — Mb) [28] in Run I. A summary of the 

mass limits is also provided in Table I. For cases (c) and (d), we simply assumed the 
three final states of bevxi, bn&Xi, an d brvx\ have the same branching ratio of 33.3%. 



TABLE I. Discovery reaches on and expected in Run II [24]. The Run I limits 

are after taking into account LEP2 limits on M~o, Mr, and M~±. In Run II, ibjfJr and 

Xi Xi 

cc^r final states will be explored by Trigger B2. 

Decay Subsequent Decay Final State of Discovery Reach in or 



(Br = 100%) Ml or Mi @20 ftT 1 (Run I) 



bi 


-&x? 


bbfir 


260 GeV/c 2 


(146 GeV/c 2 [26]) 
(116 GeV/c 2 [26]) 
(140 GeV/c 2 [28]) 
(129 GeV/c 2 [28]) 


h 




CC$T 


220 GeV/c 2 
240 GeV/c 2 


h 


— > 6Zz> £ — ► ^x? 


e+e-bfir 


h 




e + e~bfir 




h 


- bxf Xf - ^ W X? 


Ibjpr and £ + rjfir 


210 GeV/c 2 


(-) 


h 




Ibjpr 


190 GeV/c 2 


(-) 



B. Testing GMSB 

The GMSB models are generally distinguished by the presence of a nearly massless 
Goldstino (G) as the LSP. The next-to-lightest SUSY particle (NLSP) decays to its 
partner and the G. Depending on the SUSY breaking scale (y/~F) , these decays occur 
promptly (s/F < a few 100 TeV) or on a scale comparable to the size of a collider 
detector (a few 100 TeV < \[F < a few 1000 TeV). For \J~F much larger than a few 
1000 TeV, the NLSP decay takes place well outside a collider detector and are not 
directly relevant to accelerator physics. Thus we consider a systematic analysis based 
on a classification in terms of the identity of the NLSP and its decay length within the 
minimal GMSB models. The models can be specified in terms of six parameters [5]: 
N m (the number of generations of messenger fields), M m (an overall SUSY mass for 
the messengers), A (the effective visible sector SUSY breaking parameter = Fs/M m ), 
C G (the ratio of the messenger sector SUSY breaking order parameter to the intrinsic 
SUSY breaking order parameter, F/F$, controling the coupling to the Goldstino), in 
addition to tan/5 and sign(/x). The NLSP decay length scales like Cq. 

The NLSP can be bino-like Xi> higgsino-like x?, fi, or q (likely t\) in minimal 
GMSB. We here choose two scenarios, bino-like x? and fi, because the signatures define 
special detector performance at CDF and D0 other than in Section III A. A summary 
of discovery reaches for various NLSP scenarios is given in Table II. 

1. Bino-like Neutralino NLSP 

For \fF greater than a few 1000 TeV, the x? G + 7 decay takes place outside 
a collider detector. In this case, x? is essentially stable on the scale of the experiment 
and escapes as missing energy. If y/~F is less than that, however, the x? decay takes 
place within the detector. Thus, two hard photons and large fir would be observed 



TABLE II. Discovery reach (5a significance) on SUSY mass for various NLSP scenarios 
in minimal GMSB models [5]. In cr column, "p" and "d" indicate prompt and displaced 
decays of NLSP, respectively, while "11" for long-lived NLSP. 7^ indicates a displaced photon. 
Si m indicates an impact parameter. 



NLSP 


Decay Mode 


CT 


Prod. 


Key Final State(s) 


Discovery Reach 
@30 fb" 1 


Brno xi 


7 + G 


P 

d 


all 
all 


■yyftr + X 

IdjjPr or 77^ + X 


340 GeV/c 2 ( X f) 
300 GeV/c 2 (xf ) 
(cr = 50 cm) 


Higgsino Xi 


(h,Z n )+G 


P 
d 


all 
all 


(hh, hj, hZ, Z7, ZZ, 77)^r + X 
5 ip < for h -► bb, Z -► l+£- 


220 GeV/c 2 (xf") 






d 


all 


ld + X 




f 


T + G 


P 
11 


all 
all 


llljfr, PPjjfr, r h T h $r 230 GeV/c 2 (xt) 

120 GeV/c 2 (fi) 

fi(dE/dx) + ££(M« > 150 GeV/c 2 ) 420 GeV/c 2 (xf ) 

210 GeV/c 2 (f) 






11 
11 


all 
all 


n{dE/dx) + X 
/x(dE/dx+TOF) + X 


180 GeV/c 2 (f) 
210 GeV/c 2 (f) 


I co-NLSP 


e + G 


P 
d 


all 
all 


^ + dE/dx 


360 GeV/c 2 (xf) 
160 GeV/c 2 (£) 
480 GeV/c 2 (xf ) 
206 GeV/c 2 (£) 


h 


(c, bW) + G 


P 


Mi 


ccpr or £ + jets + $r 


175 GeV/c 2 (tx) 



in all final states of pair produced SUSY particles with cascade decays through the Xi 
decay. There is essentially no SM background. 

For a representative study for 77^ + X, models are chosen with N m = 1, M m /A = 
2, tan (3 = 2.5, and // > [5]. Here 77 are either prompt or displaced photons. In the 
models, the cross-section for XX (XiXi an d xtx.2) production is the largest. Thus, the 
xt is probably the best figure of merit for the discovery reach. 

For prompt photon, CDF and D0 collaborations studied the final state with differ- 
ent kinematical cuts [29], but found similar 5a discovery reach (with 30 fb _1 ) in -^^± : 

330 GeV/c 2 for CDF and 340 GeV/c 2 for D0 [5]. 

For a displaced photon, the D0 detector can reconstruct the electromagnetic (EM) 

shower develpments in the EM calorimeters to point back to the beam line to measure 

the distance of closest approach (d ca ) to the beam axis [5]. The analysis requires at 

least one photon with d ca > 5 cm in the final state of "fdjjfir- Here 7^ indicates a 

displace photon. For A = 100 TeV (M~± ~ 250 GeV/c 2 , M~ ~ 130 GeV/c 2 ), the 

Xi Xi 
probability that a photon has d ca > 5 cm is about 40% or better for cr< 100 cm. The 

discovery reach in the xt mass ranges from 310 GeV/c 2 to 280 GeV/c 2 between cr = 



and 100 cm. On the other hand, CDF will meausue the photon's arrival time on the 
EM calorimeter to distinguish from prompt photon or cosmic-ray induced photon [30]. 



2. Stau NLSP 

The models are N m = 2, M m /A = 3, tan (3 — 15, /i > with A allowed to vary [5]: 
the fi is lighter than Xi, but Mf ~ 0.5M~±. Mr — Mf is greater than a few GeV/c 2 

1 Xl 

for all points. The XiXi an d xfx° production is dominant for M~± < 350 GeV/c 2 

and the fx fx production becomes dominant for heavier xt- 

If fx is short-lived and decays in the vicinity of the production vertex, the XX 
production, followed by the cascade decays, will arise in the final states of iiijfir and 
^^33 fr- They are studied by DO, while CDF investigates the ThThftr final state from 
fx — > tG. With 30 fb _1 , both analyses have the discovery reach of 230 (120) GeV/c 2 
in M~± ~ . 

Xl (TO 

For long-lived fx, we search for events containing at least one /z-like track with a 

large dE/dx in its tracking system. The D0 analysis choose the final state including 

two leptons with M« > 50 GeV/c 2 and is sensitive for the XX production. With 

30 fb~\ the discovery reaches in M~± (M f ) are 420 (210) GeV/c 2 . 

Xi i 

The CDF detector includes a new time-of-flight (TOF) system. With a timing 
resolution of 100 ps, we could require 4<r separation (at 400 ps), which is (5^ < 2.26 (or 
p < 235 GeV/c 2 ). This sould be compared to ^7 < 0.85 using the dE/dx technique 
alone. With 30 fb _1 , the discovery reaches in Mf are 180 GeV/c 2 with dE/dx and 
210 GeV/c 2 with dE/dx plus TOF. 

It should be noted that the same technique (dE/dx, TOF) can be used [6] to test 
Anomaly-mediated SUSY Breaking (AMSB) models [31] where the xt i s long-lived 
because of a very small mass difference between xt an d Xi- 



C. Summary 



The ultimate limit for the SUSY mass reach of a hadron collider, resulting from 
the distribution function of the constituent quarks and gluons, is ~ 25% of its collider 
energy. Thus the Tevatron approaches its limit below 450 GeV/c 2 (about 50% ^/s) in 
the discovery of the SUSY particles. This would be similar to the CERN Super Proton 
Synchrotron (SppS) that approached its limit in the discovery of the W and Z bosons. 



IV. DETECTION OF B s 

We consider now the possibility of detecting the decay B s — > fi + fi~ by the CDF and 
DO detectors at the Tevatron in Run II [32]. Both detectors have been upgraded with 
excellent tracking and muon detector systems [33]. The dimuon trigger is the key to 
collect the B s — > decays. 



This process is particularly interesting for several reasons: The SM branching ratio 
is quite small, i.e., Br[B s — > /! + /i~]sm — 3.5x 10~ 9 [33]. The SUSY contribution [34-40] 
has terms that grow as tan 6 (5 and thus can become quite large for large tan (3. Finally, 
as we shall show below, the two collider detectors will be sensitive to this decay for 
tan/3> 30 in Run II. 

In order to estimate the limits on Br[B s — > /x + /x~] detection, we use the 95% 
C.L. limit on Br[B s — > fi + fJ>~] published by CDF [41]. Thus our discussion is based 
on the CDF detector, although both CDF and D0 detectors should have a similar 
p erf romance. 

In the Run I analysis, CDF observed one candidate that was consistent with B s — > 
with an estimate of 0.9 backound (BG) events in 98 pb _1 [41]. The primary Run- 
I selection variables and cut values were cr = L^MbIpj 1 > 100 /im, / = Pt^/Ipt* + 
T,p T ] > 0.75 for the muon pair, and A$ < 0.1 rad. Here, L xy is the transverse decay 
length; p 1 ^ is the transverse momentum of the dimuon system. Sp T is the scalar sum 
of the transverse momenta of all tracks, excluding the muon candidates, within a cone 
of AR = yj (Arj) 2 + (A0) 2 = 1 around the monetum vector of the muon pair. The z 
coordinate of each track along the beam line [42] must be within 5 cm of the primary 
vertex. A$ is an opening azimuthal angle between p^ and the vector pointing from the 
primary vertex to the secondary vertex (the reconstructed B- meson decay position). 
As a conservative estimate, CDF took the one event as signal to calculate 95% C.L. 
limit of signal events (Nf 5% = 5.06 events [41]) and had set a limit of Br < 2.6 x 10~ 6 . 
In the analysis, the selection efficiency (e) for signal events and the rejection power (TV) 
for background events (pass a baseline selection [41]) are estimated to be e\ = 0.45 and 
TZi = 440 by using a sample of like-sign dimuon events (5 < M w < 6 GeV/c 2 ). 

In Run II, a dimuon trigger in Ref. [33] will improve the acceptance for signal events 
by a factor of 2.8. The trigger will soon be tested using the Run Ha (2 ftT 1 ) data. This 
will allow us to modify the trigger design for the higher luminosity expected in Run 
lib (15 fb" 1 ). In this paper, we assume that the dimuon trigger can be designed by 
maintaing the acceptance for signal events. We expect to improve the acceptance for 
signal events by a factor of 2.8 [33]. If we assume the factor 2.8 to be the same for BG 
events, then we would observe 51 (386) events in 2 (15) fb _1 with the same cuts as in 
Run I. Therefore, CDF has to require a set of tighter cuts to obtain the best possible 
upper limit. 

Two types of backgrounds must be taken into account: (i) non-6 backgrounds 
comming from the primary vertex; (ii) b background events, such as the gluon-spliting 
bb events. 

One way to reduce prompt background is to require a minimum decay length L xy . 
However, two tracks can appear to form a secondary vertex if one of two tracks origi- 
nates from the primary vertex and the other has an impact parameter (S). Therefore, 
the requirement of a minimal impact parameter of individual tracks can further clean up 
the sample. It has been shown for example, in the Run-I analysis for B° — > K°*fi + fi~ 
events [43], that a tight impact parameter cut on significance for individual track 
(5 /as > 2) significantly improve the background rejection even with L xy > 100 /im. 
One has then e ~ e x x 0.43 and TZ ~ TZ± x 190. Thus a higher track impact parameter 



is neccessary to reduce the non-6 backgrounds. We would expect larger reduction with 
good efficiency even after the L xy cut. The silicon vertex detector (SVX-II) will provide 
us much better reduction for the non-6 background than Run I. 

The most severe background in Run II will be the two muons from gluon-spliting 66 
events. Since both particles are b quarks, the impact parameter does not help. Both b 
and 6 also go in the same direction, so that cut on L xy does not help either. However, 
A<3> is still usefull to remove the background events. Furthermore, in Run II, we can 
use AG in r-z view since we have z-strips in SVX-II. 

There is some room to improve the isolation cut. We can form a new isolation 
parameter by only using the tracks with large impact parameter. This new isolation 
cut will work to reject the bb rather than non-6 background. Furthermore, we can search 
for tracks with large impact parameter on the opposite side of the dimuon candidates 
to make sure that the 6 and 6 go to the opposite side. 

Therefore, CDF could improve the BG rejection by a factor of 200-400 with further 
reduction of the signal efficiency by a factor of 2-3. Based on these facts, we now 
consider two cases to evaluate Run II limits as a function of luminosity. 

In the first case (Case A), we naively assume new tighter cuts in Run II, described 
above, will gain additional BG rejection power of 450 for additonal efficiency of 0.45, 
or TZ 2 = 450 a45/e2 . This gives us 



e l l + log^/ftO/log^O) 

If we could optimize the BG rejection in Run Ha (2 fb _1 ) to be 1Z 2 ~ 5172-1 with 
e 2 ~ 0.61ei (from Eq. 1), then we would expect one BG event in 2 fb _1 . Thus, with an 
assumption of the same size of the total systematic uncertainty in Run II as in Run I, 
we can extrapolate the 95% C.L. limit to be Br < 7.7 x 10~ 8 for 2 fb -1 using Nf 5% . 

In the second case (Case B), we simply assume the Run-II background rejection 
could be improved (without loosing the signal efficiency) to keep the expected BG 
events in 2 fb" 1 at the level of Run I (i.e., 0.9 events). If we would observe one 
event in 2 fb _1 , then we could set the limits by scaling the Run-I Br limit down by 
the luminosity (2000 pb _1 /98 pb" 1 ) and the acceptance by (2.8/1.0). Thus we obtain 
Br < 4.6 x 10~ 8 . This would certainly be the optimistic scenario, but it would be a goal 
of this analysis in Run Ha. Here, the systematic uncertainty in Run II is assummed to 
be the same as in Run I. 

We repeat the same argument for different luminosities. Figure 3 shows 95% C.L. 
limits on Br[B s — > fi + fi~] at CDF in Run II as a function of integrated luminosity for 
Cases A and B. For 15 fb _1 in case A, CDF would be sensitive to Br > 1.2 x 10~ 8 and 
the combined CDF and D0 data (30 fb -1 ) would be sensitive to Br > 6.5 x 10 -9 . 

We examine first the parameter region for the mSUGRA model that would be 
accessible to CDF or D0 at Run II with 15 fb _1 of data. Figure 4(left) shows the 
Br[B s — > fi + fi~] as a function of mi/ 2 for A = 0, m = 300 GeV. One sees that with a 
sensitivity of Br[B s — > fi + > 1.2 x 10~ 8 for 15 fb _1 , the Tevatron Run II can probe 
the B s — > fi + fi~ decay for tan/3 > 30. Further, a search for this decay would sample 
much higher regions of mi/ 2 than a direct search at Run II for SUSY particles which 
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FIG. 3. Illustrated 95% C.L. limits on the branching ratio for B s — ► /U + /U _ at CDF in Run 
II as a function of integrated luminosity [32]. Solid (Case A) and dashed (Case B) curves are 
based on different assumptions on the signal selection efficiency and the background rejection 
power. See the text for details. 



is restricted to mi/2 < 250 GeV [13]. As m increases, the branching ratio goes down. 
However, this dependence becomes less significant for large m^, where m as large as 
800 GeV can be sampled for large m^. 

In Figure 4(right) the contours of Br[B s — > [i + fi~] are plotted in the m -mi/ 2 plane 
for tan/5 = 50, A = 0. Also see Figure 1. We combine now this result with the 
other experimental constraints. Thus the shaded region to the left is ruled out by 
the b — > S7 constraint, and the shaded region on the right hand side is disallowed if 
fl suGRA > 1 1 x io~ 10 . The narrow shaded band in the middle is allowed by the dark 
matter constraint. We note that independent of whether the astronomically observed 
dark matter is SUSY in origin, the dark matter allowed region for mSUGRA cannot 
significantly deviate from this shaded region, for below the narrow shaded band, the 
stau would be lighter than the neutralino (leading to charged dark matter), while above 
the band, mSUGRA would predict more neutralino dark matter than is observed. 

Using our estimate that Br > 1.2 x 10~ 8 can be observed with 15 fb -1 , we see that 
almost the entire parameter space allowed by the < 11 x 10~ 10 constraint can be 
probed in Run II for tan/3 = 50. Note that an observed Br[B s — > /i + /i~] > 7 x 10~ 8 , 
possible with only 2 fb _1 (see Figure 3), would be sufficient to rule out the mSUGRA 
model for tan/3 < 50. In Figure 4(right) we also show the expected dark matter 
detector cross section for Milky Way dark matter (the short solid lines). For tan/3 = 
40 and Aq = 0, we see about half the parameter space can be scanned by the CDF 
detector (and the whole parameter space with 30 fb" 1 if a similar analysis holds for the 
D0 detector). We note, further, that if A = 0, a simultaneous measurement of both 
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FIG. 4. Branching ratio for B s — > in mSUGRA models [32]. On the left plot, 

Br[B s — > n + fJ.~] is shown as a function of m 1 / 2 for various tan/3 values, where dashed and 
dash-dotted lines are to indicate the models are excluded via constraints on Br[b — > 57] 
and rrif > m~o, respectively. On the right plot, three dashed lines from left to right in the 

mQ-mi/2 plane indicate Br = 7 x 10 -8 , 2 x 1CP 8 , 1 x 10~ 8 , respectively, for tan/3 = 50. The 
three short solid lines indicate the a^o_ p values (from left: 0.05 xl0~ 6 pb, 0.004 xlO -6 pb, 
0.002 xl0~ 6 pb). The vertical dotted lines label Higgs masses. 



Br[B s — > and a M would essentially determine the mSUGRA parameters, as the 

mo allowed region at fixed mi/2 is very narrow due to the dark matter constraint. The 
effect of varying A is also discussed in Ref. [32]. 

In Figures 4(right), we have also drawn lines for various light Higgs masses (ver- 
tical dotted lines). A measurement of B s — > fi + fi~, and would then effectively 
determine the parameters of mSUGRA for > by requiring that they intersect with 
the dark matter allowed band at a point. (If no choice of parameters allowed this, 
mSUGRA would be ruled out.) The Tevatron Run II should be able to either rule out 
a Higgs mass or give evidence for its existence at the 3a level over the entire allowed 
mass range of SUSY light Higgs masses. Alternatively, the LHGs determination of mh 
or the gluino mass (to determine m^j) would fix the parameters of mSUGRA. 

Thus if a^i increases, the bound moves downward, encroaching further on the allowed 
part of the parameter space, and a value of a M > 50 x 1CU 10 would eliminate the 
mSUGRA model [44]. However, if a M decreases significantly (but is still positive), 
the mSUGRA model would predict a heavy SUSY particle spectra closer to the TeV 
region, having significant effects on accelerator and dark matter detection physics. An 
accurate determination of a M corresponds to a line from upper left to lower right (or 
more precisely a band when errors are included) running parallel to the < 11 x 10~ 10 



boundary, and cutting through the allowed dark matter band which runs from lower 
left to upper right. Thus, these two experiments are complementary for determining 
the mSUGRA parameters. 

It is also pointed out that the B s — > decay is very powerful in testing the 

SUSY breaking mechanism [45]. If the B s — > decay is observed with Br > 1CT 8 , 

the Tevatron could exclude (i) GMSB models (with tan/3 < 40 and lower N m ) and 
(ii) minimal AMSB models. 

Interestingly, unlike the R parity conserving models, the SUSY contribution to 
B s — > can occur at the tree level in the R parity violation models. The 

Br[B s — » can be large for both small and large tan/3. We, for example, consider 

lepton violating terms [46], and so we set A" to zero (to prevent rapid proton decay). 
In Ref. [32], the branching ratio of (feasible with 2 fb _1 ) could be sensitive upto 
mi/a < 1000 (750) GeV/c 2 for m = 300 (500) GeV/c 2 . 

V. CONCLUSION 

We have investigated the prospects for studying SUSY parameter space in Run II 
of the Fermilab Tevatron. Two complementary obeservables (especially in SUGRA 
models) are of great interest [15]: trilepton final state as production of charginos and 
neutralinos for lower tan (5 values and the rare decay mode B s — > for large tan (3 

values. 

Especially, the B s — > ^ decay is an important process for SUSY searches, as the 
Standard Model prediction of the branching ratio is quite small (3.5 x 10~ 9 ), and the 
SUSY contribution increases for large tan/3 as tan 6 /3. We find that a Br > 1.2 x 10~ 8 
for 15 fb" 1 can be probed by each collider detector in Run II. (This is nearly a factor 
of 100 improvement over the Run I bound.) For the mSUGRA model, the above 
sensitivity implies that Run II could probe a region of parameter space for tan/5 > 30, 
a region which could not be probed by a direct search at Run II. A large branching 
ratio, i.e., > 7 x 10~ 8 (14 x 10~ 8 ) would be sufficient to eliminate the mSUGRA model 
for tan/? < 50(55) if one combined the expectations for B s — > fi + fi~ for mSUGRA with 
other experimental bounds on the parameter space. 
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